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Abstract 



pi ' X-ray and near-infrared {J-H-Ks) observations of the Galactic black hole binary GX 339-4 in the 

low/hard state were performed with Suzaku and IRSF in 2009 March. The spectrum in the 0.5-300 keV 
band is dominated by thermal Comptonization of multicolor disk photons, with a small contribution from a 
direct disk component, indicating that the inner disk is almost fully covered by hot corona with an electron 
^ ■ temperature of wl75 keV. The Comptonizing corona has at least two optical depths, r 1,0.4. Analysis 

of the iron-K line profile yields an inner disk radius of (13.3lg'o)i?g (i?g represents the gravitational radius 
GM/c^), with the best-fit inclination angle of w 50°. This radius is consistent with that estimated from 
^ I the continuum fit by assuming the conservation of photon numbers in Comptonization. Our results suggest 

. that the standard disk of GX 339-4 is likely truncated before reaching the innermost stable circular orbit 

OO ' (for a non rotating black hole) in the low/hard state at ^1% of the Eddington luminosity. The one-day 

•/^ , averaged near- infrared light curves are found to be correlated with hard X-ray flux with Fks F^^^. 

' The flatter near infrared vF^, spectrum than the radio one suggests that the optically thin synchrotron 

If^ I radiation from the compact jets dominates the near-infrared flux. Based on a simple analysis, we estimate 

' the magnetic fleld and size of the jet base to be 5 x 10* G and 6 x 10^ cm, respectively. The synchrotron 

self Compton component is estimated to be approximately 0.4% of the total X-ray flux. 

Key words: accretion, accretion disks — black hole physics — infrared: stars — stars: individual (GX 
339-4) — X-rays: binaries 

1. Introduction suggest that the disk is extending down to the innermost 

stable circular orbit (ISCO). 

Black hole binaries (BHBs), a class of the brightest In contrast, at lower accretion rates, the hard X-ray 
Galactic X-ray sources, are ideal objects to study the ac- (> 10 keV) contribution to the entire flux becomes more 
cretion disk physics around a black hole. They exhibit substantial. This state, called the "low/hard" state, shows 
distinct "states" with different luminosity and spectral X-ray spectra approximated by a power law of photon in- 
shape, which are mainly determined by the mass accre- dex 1.5 — 1.9 with a high energy cut-off at ^ 100 keV. 
tion rate (see e.g.. Done et al. 2007 for a recent review). This emission is interpreted to be produced by unsat- 
When the accretion rate is high (but below several 10% of urated Compton up-scattering by hot electrons of seed 
the Eddington limit), the spectra of BHBs are dominated photons emitted by the disk. The physical properties and 
by a soft X-ray component below 10 keV. This state is geometry of the corona in relation to the disk, however, 
called the "high/soft" state, and its dominant soft com- remain to be studied. While there is a consensus that 
ponent is described by a Multi-Color Disk (MCD) model, an optically thick disk is truncated before reaching the 
formed by a superposition of blackbody radiation at dif- ISCO at low mass accretion rates (e.g., Esin et al. 1997), 
fercnt temperatures emitted from various radii of the disk discussion at what Eddington fraction the disk becomes 
(Mitsuda et al. 1984). In this state, many observations truncated is still controversial (see below). 
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Radio to optical observations of BHBs have provided ev- 
idence for the presence of stable jets in the low/hard state, 
while the jet is quenched in the high/soft state (Fender ct 
al. 1999a). By using radio interferometers, the central 
compact core and jet-like knots are resolved in several 
BHBs in the low/hard state (Fender 2001, and references 
therein). In the radio region, the spectral energy distribu- 
tion (SED) of BHBs exhibits an almost flat profile, which 
is considered to be optically thick synchrotron radiation 
from the jets, like those observed from active galactic nu- 
clei (Blandford & Konigl 1979). This flat spectrum is 
thought to extend to the sub-mm or infrared wavelengths, 
above which the synchrotron emission becomes optically 
thin and eventually declines as a power law with an index 
> 1 (SED units). Thus, observations at the wavelengths 
where the synchrotron luminosity from the jets peaks give 
key information to understand the energetics of the jet. 

GX 339-4 is a transient Galactic BHB, discovered in the 
early 1970s (Markert et al. 1973). Ahhough this object 
was observed for many years at various wavelengths, its 
binary-system parameters arc still uncertain (see Section 
2). The accretion disk geometry in the low/hard state 
has been studied from X-ray spectra at different lumi- 
nosities. Tomsick et al. (2009) show that, at ^ 0.001 of 
the Eddingtion luminosity (iEdd), the standard disk is 
truncated at > 35i?g (i?g = GM/(? is the gravitational 
radius, where G, M, and c are the gravitational constant, 
black hole mass, and light velocity, respectively.) Miller 
et al. 2006 argue that the disk extends to ~ 4i?g in the 
bright low/hard state at ~ 0.06 of the Eddington lumi- 
nosity, based on the detection of a broad iron-K emission 
line. Similar results are reported by Tomsick et al. (2008) 
at 0.023iEdd and O.OOSLEdd- Done & Diaz Trigo (2010), 
however, re-analyze the same data used by Miller et al. 
(2006), and find that the spectra are strongly affected by 
pile- up (multiple photon events in CCD). Extracting the 
simultaneous data that are free from pile-up, they obtain 
a narrower iron-K line profile that rather favors a trun- 
cated disk with an innermost radius of > 6i?G- Thus, the 
disk evolution of GX 339-4 in the low/hard state as a 
function of luminosity is not firmly established yet. GX 
339-4 is also a persistent radio source in the low/hard 
state, whose spectrum is flat or slightly increases accord- 
ing to the frequency with the spectral index « -1-0.1-0.2 
(Corbel et al. 2000). In spite of the importance of the 
near-infrared data, simultaneous multi-band observations 
in the J-H-K bands and X-rays have been very limited 
so far. 

To establish the inner disk structure of BHBs in the 
low/hard state, we observed GX 339-4 using the X-ray 
satellite Suzaku in 2009 March, when it became active and 
stayed in the low/hard state. Quasi-simultaneous near- 
infrared observations were performed with the Infrared 
Survey Facility (IRSF) 1.4 m telescope to reveal the jet 
activity. In this paper, we first summarize the constraints 
on the system parameters of GX 339-4 based on previous 
observations in Section 2. The Suzaku observations and 
data reduction are described in Section 3, and analysis 
of the X-ray data is given in Section 4. In Section 5, we 



present details of the IRSF observations and data analysis. 
We discuss our X-ray and near-IR results by comparing 
with previous observations of GX 339-4 and other BHBs 
in Section 6. The conclusions are summarized in Section 
7. 

2. Summary of Constraints on Inclination and 
Distance of GX 339-4 

The system parameters of GX 339-4 including the dis- 
tance, inclination angle, black hole mass (hereafter A/bh), 
and the mass of companion star {Mq) are still unknown. 
Figure 1 summarizes the constraints derived from previous 
observations in the inclination (or black hole mass) ver- 
sus distance plane. The limit that i is less than ^ 60° is 
imposed by the fact that GX 339-4 exhibits no eclipse 
(Cowley et al. 2002). The upper and lower limits of 
the distance D are estimated in several different ways. 
Zdziarski et al. (1998) found that D > 3 kpc by investi- 
gating the interstellar reddening, while Hynes et al. (2004) 
estimated 6 kpc < D < lb kpc using the structure of 
NaD line obtained through high-resolution spectroscopy. 
Zdziarski et al. (2004) put the hmit 7 kpc < I? < 9 kpc, 
using the apparent radius of the secondary low-mass star 
based on optical and infrared observations. In Figure 1, 
we plot the Hynes et al. (2004) and Zdziarski et al. (1998) 
results for the upper and lower limits of the distance, re- 
spectively. 

The dependence of TI/bh on the inclination angle is ob- 
tained from the mass function, which is expressed as 
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where P, G, and Kq are the orbital period, the gravita- 
tional constant, and the radial velocity semi- amplitude of 
the companion star, respectively. Hynes et al. (2003) es- 
timated /(Mbh) of GX 339-4 to be 5.8 Mq (2.0 Mq at 
a 95% confidence lower limit) by observing the Doppler 
shifts of the N III and He II emission lines from the irra- 
diated companion star. Ignoring the Mq/Mbk term in 
Equation 1, whose value is less than 0.08 according to 
Zdziarski et al. (2004), we calculate the black hole mass 
corresponding to each inclination angle, which is indicated 
in the upper axis. 

For a given /(Mbh) value, we can also constrain the re- 
lation between the distance and inclination angle, assum- 
ing that the fraction of the Eddington luminosity in the 
very high state should be in the range of 0.1 ^L/L-^dd ^ 1- 
Here we refer to the 0.2-100 keV luminosity in the very 
high state obtained with Suzaku in 2007 (Yamada et al. 
2009), 3.8 X 10^^{D/8 kpcf erg s'^. 

Using the proper motion of the relativistic jets emitted 
from the system, we are also able to connect the distance 
and inclination angle (assumed to be perpendicular to the 
jet axis) for an assumed intrinsic speed of the jet. In 2002, 
bright radio knots are observed from GX 339-4 at the 
state transition and the head of one knot moved 12 arcsec 
in approximately 300 days (Gallo ct al. 2004). From this 
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Fig. 1. Relation of the inclination angle, black hole mass, and distance of GX 339-4. The abscissa represents the inclination angle 
(bottom) and black hole mass (top). The mass function is assumed to be 5.8 Mq in the left panel and 2.0 Mq in the right panel. 
The orange horizontal line shows the upper limit of distance (15 kpc) derived from the velocity dispersions (Hynes et al. 2004), and 
the purple line corresponds the lower limit (3 kpc) determined with the extinction study (Zdziarski et al. 1998). The vertical line at 
i = 60° shows the upper limit of the orbital inclination angle by the absence of eclipse (Cowley et al. 2002). The relation between 
the inclination angle and distance, calculated from the 0.5—200 keV flux of GX 339—4 in very high state (Yamada et al. 2009) is 
presented by assuming L/LEdd = 1 (pink, upper curve) and 0.1 (dark blue, lower curve), with the area for 0.1 < L/LEdd < 1 filled in 
gray. The red and green dotted lines show the relation between the inclination and distance derived from the jet motion observed 
by Gallo et al. (2004), where the intrinsic velocity is assumed to be 0.92c and 0.98c, respectively. 



apparent velocity, we derive the relation between the incli- 
nation angle and distance, assuming the intrinsic speed of 
the jet to be 0.98c (Fender et al. 1999b) or 0.92c (Mirabel 
& Rodriguez 1994), a typical value observed from the mi- 
croquasar GRS 1915-f 105. 

Also, recent work sets the lower limit of the companion- 
star mass, Mq > 0.166A/q, by considering the mass trans- 
fer process of low mass X-ray binaries (Munoz-Darias et 
al. 2008). This minimum mass makes a tight constraint on 
the binary inclination, i > 45, by assuming that the black 
hole mass is less than 15M0, which is the largest mass ever 
known among Galactic BHBs (discussed in Kolehmainen 
& Done 2010). This is quite consistent with the con- 
straints summarized in Figure 1 for both cases of the two 
/(Mbh) values. 

3. X-Ray Observations and Data Reduction 

3.1. Observations 

We performed three sequential ToO (Target of 
Opportunity) observations of GX 339-4 with Suzaku 
(Mitsuda et al. 2007) on 2009 March 18, 25-26, and 30- 
31 (hereafter Epoch-1, -2, and -3, respectively), each for a 
net exposure of ~40 ksec. The ToO observations were 
triggered after the increase of the X-ray flux detected 
by RXTE/PCA scanning observations in 2009 February 
(Markwardt et al. 2009) and by Swift/BAT monitoring in 
March (Still et al. 2009), based on our approved Suzaku 
AO-3 program. 

In Figure 2 we show the Suzaku hard-state observations 
in the context of the large data archive for GX339-4 pro- 
duced with the RXTE PCA Instrument. We extracted 
the normalized PCA count rates in four energy bands, as 
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Fig. 2. Hardness-intensity diagram (upper) and color-color 
diagram (lower) of GX 339-4, obtained with RXTE/PCA 
during 1996-2010. The four orbits during our Suzaku cam- 
paign is highlighted in red (filled triangle), one point on MJD 
54914 and three on MJD 54916. 
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Table 1. Summary of the Suzaku observations 

Label Observation Date (UT) Observation ID Net Exposure (ksec) 

Start End XIS-0 and -3* XIS-l^ HXD 

Epoch-l 2009 March 18 01:53:13 2009 March 18 23:04:19 403011010 21.4 3.9 35.4 

Epoch-2 2009 March 25 08:14:00 2009 March 26 06:30:24 403011020 19.4 3.6 34.8 

Epoch-3 2009 March 30 11:27:27 2009 March 31 12:31:10 403011030 19J 3.6 35.1 

* XIS-0 and -3 were operated in the 1/4 window mode with 0.5 see burst option, 
t XIS-1 was in the full window mode with 0.5 sec burst option. 



described by Lin et al. (2007), to construct the hardness- 
intensity diagram (HID; top panel) and color-color dia- 
gram (bottom panel). Each of the 1159 plotted points 
represents an interval of continuous exposure, with an av- 
erage exposure time of 2.0 ks. These data have been fil- 
tered to eliminate exposures obtained when the source flux 
is fainter than 12.5 count sec""'^ PCU~^ (or 5 mCrab at 
2-30 keV), so that the statistical error bars remain similar 
to or smaller than the size of the plot symbols. 

Four of these RXTE exposures overlap with the Suzaku 
observations (one on MJD 54914 and three on 54916), and 
the results for those data are plotted as red triangles. The 
HID in Figure 2 extends the results shown by Remillard 
& McClintock (2006), Figure 9d to cover RXTE observa- 
tions through 2010 July 8. In the earlier work, the HID 
points are symbol-coded to show state assignments based 
on quantitative definitions, and observations in the hard 
state completely surround the regions occupied by the red 
points (i.e., the Suzaku-associated exposures) shown in 
Figure 2. 

The Suzaku observations of GX 339-4 were conducted 
by using the X-ray Imaging Spectrometers (XISs) and 
the Hard X-ray Detector (HXD) at the HXD nominal 
position. The XISs arc composed of three frontside- 
illuminated (FI) CCDs (XIS-0, -2 and -3) and a backside- 
illuminated (BI) CCD (XIS-1) covering the energy range 
of 0.2-12 keV. The BLCCD has larger effective area than 
FI-CCDs below «1.5 keV and is more sensitive to low en- 
ergy photons. The HXD consists of silicon PIN diodes and 
GSO crystal scintillators, and is sensitive to hard X-ray 
photons in 10-70 keV and 40-600 keV bands, respectively. 
In our observations, the editing mode of the FI-CCDs was 
2x2 and that of the BI-CCD was 3x3. To reduce photon 
pile-up, XIS-0 and XIS-3 were operated with the 1/4 win- 
dow and 0.5 sec burst option, and XIS-1 was used with 
the full window and 0.5 sec burst option. XIS-2 was not 
usable then because of the trouble of the instrument on 
2006 November 9. The summary of the observations is 
given in Table 1. 

3.2. Data Reduction 

We utilized data products of both instruments reduced 
through Suzaku pipeline processing version 2.3.12.25. To 
extract light curves and spectra, the data were analyzed 
with HEAsoft version 6.6.3 and the calibration data files 
(CALDB) released on 2009 October 7. Because the 
Suzaku attitude was still unstable shortly after the ma- 
neuver to the target, the data acquired in the first 2 ksec 
in Epoch-l were discarded to avoid difficulties of analysis. 



We analyzed the XIS cleaned event files in the "sky coor- 
dinates" system, where the small attitude variation of the 
spacecraft during the orbit is corrected. We extracted the 
counts of the FTCCDs (XIS-0 and -3) and the BI-CCD 
(XIS-1) from circular regions centered on the target with 
the radii of 1.9' and 3.7', respectively. The regions enclose 
^90% of the available photons in each field of view, whose 
area is a factor of 4 smaller in FI-CCDs than in BI-CCD 
because of the window option. The background of the XIS 
were taken from a source-free area with the same radii as 
the target. For the HXD, we used the modeled back- 
ground files produced by the Suzaku HXD team^. The 
PIN background, composed of the non X-ray background 
(NXB) and the cosmic X-ray background (CXB), was sub- 
tracted from the source spectra. Although there is addi- 
tional background contribution from the Galactic ridge 
emission, the effect is estimated to be negligible (w 0.5% 
of the target flux) in our PIN data. For the background of 
GSO, only the NXB component was considered because 
the contribution of the CXB to the total background rate 
is less than 0.1% ^. 

The light curves of GX 339-4 in four energy bands, 0.7- 
2 keV, 2-10 keV, 14-50 keV, and 50-200 kcV are plotted 
in Figure 3 for the three epochs. The PIN and GSO light 
curves were corrected for dead time, and the NXB was 
subtracted. Background subtraction was also performed 
for the XIS data. As presented in Figure 3, the count 
rates of the target in these energy bands exhibit a sim- 
ilar behavior to one another, suggesting that the overall 
spectral shape is mostly constant. 

Figure 4 shows the time-averaged spectra of the XIS 
and the HXD in the three observations. To improve 
statistics, XIS-0 and -3 were combined each other. We 
created the response matrix files and auxiliary response 
files of the XIS with the FTOOLS xisrmfgen and 
xissimarf gen, respectively. As the response files of 
the HXD, we used aejixd_pinhxnome5_20080716 . rsp for 
PIN, and ae_lixd_gsohxnomJ20080129 . rsp and aeJixd_ 
gsohxnom_crab_20070502 . arf ^ for GSO. To include pos- 
sible uncertainties associated with the calibration of the 
instruments, 1% systematic error is introduced for each 
spectral bin of the XIS and the HXD. 

To examine possible effects of pile-up in the XIS data, 
which should be more significant in the core region of the 
point spread function of the X-ray telescope, we also make 

^ http: / /www. astro . isas .jaxa.jp/ suzaku/analysis/hxd/pinnxb , 
gsonxb 

^ http:/ /www. astro. isas. ac.jp/suzaku/analysis/hxd/gsonxb/ 
http:/ /www. astro. isas.jaxa.jp/suzaku/analysis/hxd/gsoarf/ 
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Fig. 3. X-ray light curves of GX 339-4 binned in 5760 sec. The ordinatcs represent the count rate (counts sec ^) in the energy 
range of (a) XIS-0 0.7-2 keV, (b) XIS-0 2-10 keV, (c) PIN 14-50 keV, and (d) GSO 50-200 keV. The left, middle, and right panels 
present the results of Epoch-1, -2, and -3, respectively. The background subtraction (for the XIS and HXD) and dead time correction 
(for the HXD) are applied. 




Epoch 1+2+3, XIS0+XIS3 



Fig. 4. Time-averaged Suzaku spectra obtained in Epoch-1 
(black, open square), -2 (red, open triangle) and -3 (blue, 
filled circle). For clarity the XIS-1 data are not plotted. 



Energy (keV) 

Fig. 6. Ratios of the XlS-O-l-XlS-3 spectra (sum of Epochs 1, 
2, and 3) extracted from the whole circular region (r < 1.9', 
black) and annulus of r = 0.57' — 1.9' (blue, open square), 
divided by that from the r = 0.95' — 1.9' annulus. They are 
normalized to be unity at 5.0 keV. 



00f 



Fig. 5. XIS-0 image of Epoch-1 in the sky coordinates. 
The circles (red) correspond to radii of 1.9', 0.95', and 0.57' 
(dashed), which define the annulus regions used to check the 
pile-up effects. 



the spectra of XIS-0 and XIS-3 in two annuli whose inner 
radii are set to be 30% (0.57') and 50% (0.95') of the outer 
radius (1.9'). The regions are ihustrated in Figure 5 for 
XIS-0. For each annulus region, both XIS-0 and XIS-3 
spectra in the three epochs are combined into a single 
one. Figure 6 plots the ratios of the spectra extracted 
from the whole circular region (r < 1.9', black) and the 
r = 0.57' — 1.9' annulus (blue), divided by that from the 
r = 0.95' — 1.9' annulus, which can be regarded as free from 
pile-up. As noticed from the figure, the spectral ratio for 
the whole circular region shows a significant excess above 
~8 kcV attributable to pile-up. Thus, while we adopt the 
r < 1.9' spectra to achieve the best statistics, the energy 
range above 8 keV is not utilized in our spectral analysis. 
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The difference from the r = 0.95' — 1.9' spectrum is very 
small below this energy band, and we confirm that our 
results on the innermost disk radius derived from the iron- 
K profile (Section 4.3) are not affected by pile-up. We 
find that its estimate from the continuum fit (Section 4.2) 
could be affected by ~20% in Epoch 1 and by < 10% in 
Epochs 2 and 3, but this uncertainty does not change our 
conclusion on the black hole mass presented in Section 4.3. 

4. Analysis of the Time- Averaged Suzaku Spectra 

In this section, we present the analysis of the Suzaku 
spectra extracted above. The spectral fit is performed by 
using XSPEC version 11.3.2ag. In all spectral models, 
we assume solar abundances given by Anders & Grevesse 
(1989). We employ the wabs model for Galactic interstel- 
lar absorption, whose column density is set to be a free 
parameter. The quoted errors refer to 90% confidence lev- 
els for a single parameter. 

4.I. Broad Band Fitting 

To best constrain the origin of the X-ray emission, the 
broad-band spectra of the XIS and HXD are simultane- 
ously analyzed. We use the data in the energy range of 
0.5-8.0 keV for XIS-1, 1.0-8.0 keV for XIS-0 plus -3, 15-50 
keV for PIN, and 50-310 keV for GSO. To avoid the cal- 
ibration uncertainties around the instrumental Si-K and 
Au-M edges in the responses, the spectra in the 1.7-1.9 
keV and 2.1-2.3 keV energy band are excluded in the fit. 
The cross normalization between HXD and XlS-O-fXIS- 
3 spectra is fixed at 1.18 based on the calibration result 
using the Crab Nebula data ^, while that between XIS-1 
and XIS-04-XIS-3 is set free. 

As shown in Figure 4, the spectra of the three observa- 
tions exhibit similar shape to one another, indicating that 
the physical state of the source remain the same in the 
three epochs. We find that the overall continuum in the 
2-100 keV band can be roughly approximated by a power 
law of a photon index of ~1.5 in each epoch. This con- 
firms that GX 339-4 was the low/hard state throughout 
our observations. 

Following the previous studies of BHBs in the low/hard 
state with Suzaku, we adopt a model from a standard 
accretion disk and its thermal Comptonization by a hot 
corona as the basic description of the broad band contin- 
uum. The disk emission is assumed to be the multi-color 
blackbody radiation modeled by diskbb (Mitsuda et al. 
1984). For the Comptonization component, we employ 
the compPS model (Poutanen & Svensson 1996), where 
the Comptonized spectrum is computed from the electron 
temperature, the optical depth of scattering, and the en- 
ergy distribution of incident photons. We consider only 
thermal electrons by setting gmin= —1. We assume that 
the Comptonization cloud is spherical (i.e., geom=4) and 
that seed photons are entirely originated in the standard 
disk (Tbb = -Tin in the XSPEC terminology). 

The compPS code is able to take into account a reflec- 
tion component from the disk as calculated by Magdziarz 
& Zdziarski (1995). Its parameters are the solid angle 



of the reflector visible from the source (fi/27r), ioniza- 
tion parameter (^), inclination angle (i), and tempera- 
ture, which we fix at 10^ K. This component is further 
relativistically blurred with the diskline profile (Fabian 
et al. 1989), by assuming a radial dependence of emis- 
sivity as oc between the inner and outer disk radius 
(i?in and i?out)- We fix i?in = lORg, i?out = 10000i?g, 
i = 50°, and /3 = —2.3, based on the results from the anal- 
ysis of the XIS spectra in the 3-8 keV band (Section 4.3). 
Although Compton broadening is not taken into account 
in the reflection model, the effect is negligible in our case 
because the ionization parameter is found to be very low 
(^ < 2, Section 4.1) and estimated temperature does not 
exceed '--^ 10^ K as assumed. The iron K fluorescent line 
at ~ 6.4 keV is simply fitted with a Gaussian at this 
stage, whose detailed modeling does not affect the anal- 
ysis of the broad band continuum. Combining all the 
components described above, the model is expressed as 
wabs X (compPS + diskbb + gauss) on the XSPEC termi- 
nology. 

As the first trial, we perform a spectral fit without re- 
flection components. The result is far from acceptable, 
yielding x^/zy = 1219/950 (Epoch-1), 1160/907 (Epoch- 
2), and 1195/873 (Epoch-3), with residuals shown in the 
second panels of Figure 7. As noticed, there remains an 
absorption edge-like feature above ~ 7 keV in the XIS 
spectra and a hump structure at ^ 30 keV in the PIN 
data. They are strong evidence for the reflection compo- 
nent, as reported in previous broad band observations of 
GX 339-4 in the same state (e.g., Ueda et al. 1994). 

We find that inclusion of the reflection component in 
the model greatly improve the fit, although the fit good- 
ness is not yet satisfactory, ^lv= 1161/948, 1083/905, 
and 1136/871 for Epoch-1, -2, and -3, respectively. There 
is still a discrepancy between the data and model in the 
15-50 keV region, as shown in the third panels of Figure 
7. We confirm that this is not caused by any systematic 
uncertainties in the response or background of PIN, based 
on the calibration using the Crab Nebula and comparison 
of the background model with the dark-Earth occulta- 
tion data in our observations. Similar residuals are also 
found from other BHBs in the low/hard state when fit- 
ted with the same model (GRO J1655-40 by Takahashi 
et al. 2008 and Cyg X-1 by Makishima et al. 2008). A 
possible reason for the discrepancy is that the corona can- 
not be described by a constant temperature and optical 
depth, and therefore we add another compPS component. 
This "double compPS" model was successfully applied to 
the low/hard state spectra of BHBs, such as Cyg X-1 
(Makishima et al. 2008; Gierlihski et al. 1997) and GRO 
J1655-40 (Takahashi et al. 2008). We tie all the param- 
eters of the additional compPS component to those of the 
other compPS component except for the optical depth r 
and normalization. 

The best-fit model and residuals of fit using the dou- 
ble compPS model are displayed in the top and bottom 
panels of Figure 7, respectively. In Figure 8, the best- 
fit model spectrum in the vF^ form (where F^ is the en- 
ergy flux) is plotted, with contribution of each component. 
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Fig. 7. Simultaneous fit to the XIS-1 (red, open circle), XIS-O+XIS-3 (black:), PIN (green), and GSO (blue) spectra with the double 
compPS model in Epoch-1, -2, and -3. The second, third, and bottom panels show the residuals of the fit with the single compPS 
model without reflection, that with reflection, and the double compPS model with reflection, respectively. 
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Fig. 8. The best-flt double compPS model of each epoch in the vF^ form, corrected for the interstellar absorption. Each component 
is separately plotted. 
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Table 2. Parameters of the Double Compton Model* for Epoch-1, -2, and -3 



Component 


Parameter 


Epoch-1 


Epoch-2 


Epoch-3 


wabs 


iVfj (i^'J cm ) 


K Q _i_ n no 
O.oO ± U.Uo 


C- r9 + U.Ub 

O-"J^-0.03 


0.4D ± U.Uo 


diskbb 


Jin (keV) 


0.ZZ3 ± O.L)U3 


n 990+0.002 

0.220_o.o2 


r, 99r+U.0U2 
U.ZZO_o.o03 




normalization^ 


-'ion 
<lzO 


<zoo 


9i-n+120 




photon flux' 


i.oy X iU 


i.Ol X iU 


4.Uy X lU 


compPS^ 

(harder component) 


Jo (KeVj 
(y-parameter'^ 

T 


1 7R+4 

176_6 

1 r;c+0.04 
^•3^-0.01 
1 -19+0.07 
-■^■^^-0.03 


1.2±0.4 
0.9±0.3 


177_g 

1 24+°-^^^ 

^■^^-0.01/ 
^■^^-0.03 




photon flux" 


0.283 


0.679 


0.547 


compPS''' 


(y-parameter'^ 


0.63 ±0.01 


0.40 ±0.02 


o.siiH^) 


(softer component) 


T 


r\ Ac^ J- r\ no 
U.4d ± U.Uz 


n qn+0.02 
0.30_o.oi 


U.o / ± U.Uz 




photon flux" 


1.10 


0.889 


0.798 


renection+ 
Gaussian 


(" / — 1^ 
^ (erg cm s ) 

Scon (kcV) 

EW (eV) 
a (keV) 


0.474+0470 
0.46 ±0.01 

"•^-■^-0.06 
^■'J°-0.03 


< 1.71 

0.48t°:°? 
6.44 ±0.06 

4o 11 

0-38_o;o6 


0.481Io.477 

U-^'J-0.01 

"■^"-0.03 

0.181^;°^ 




(km)** 


117±2 


128lf 


115 ±2 




xVd.o.f. 


996/946 


873/903 


940/869 


Unabsorbed flux 


0.5-2 keV (erg cm"^ s"^) 


3.5 X 10-^" 


3.8 X 10-"' 


3.5 X 10-^" 




2-10 keV (erg cm"^ g-i) 


6.4 X 10^1" 


6.6 X 10^1° 


6.2 X 10^1° 




10-300 keV (erg cm^^ g-^) 


4.4 X lO-'^ 


4.3 X 10-^ 


4.0 X 10-9 



Notes, 'i'he errors are 90% confidence level for a single parameter. 

* The model is expressed as wabs X (compPS + compPS + diskbb + gauss). 

t We set gmin = — i (thermal electrons only) and goom = 4 (spherical geometry). The seed photons for 
Comptonization is the MCD emission from the disk (i.e., T^b in compPS is linked as —Tin in diskbb). 
t We assume the solar abundances for iron and heavy elements. The temperature of the reflector is fixed at fO® K. 
The reflection parameters of the two compPS components are linked each other, with JJin, JJout, /3, and i fixed 
at the best-fit values from the iron-K line fit summarized in Table 3. 

^ ((?'in/km)/(D/iOkpc))^ COS0, where and 9 are the inner disk radius and the disk inclination, respectively. 

II Unabsorbed photon fiux in units of photons cm"'^ s~^ in the 0.01-fOO keV band based on the best-fit parameters. 

# Compton y-parameter calculated as 4TkTe/{meC^). 

•* Estimated innermost disk radius derived from the fiux of direct MCD and Comptonization components, 
assuming a distance of 8 kpc and an inclination of 50° . Those have been corrected for the inner boundary condition 
and color /effective temperature ratio (see Section 4.2). 



corrected for the interstellar absorption. The best-fit val- 
ues are summarized in Table 2. The fit is significantly 
improved, xV^^ = 996/946, 873/903, 940/869, compared 
with the single compPS model, yielding F-test probabilities 
of 4 X 10-32, 7 X 10-43, and 1 x 10-^6 for Epoch-1, -2, and - 
3, respectively. We find that the contribution of the direct 
MCD component (diskbb) is very small (< 2% in flux at 
0.5 keV), and only its upper limit is obtained in Epochs 1 
and 2. The absorption column density is determined as 
A^H ~ 5.4 x 10^^ cm^^^ which is consistent with previous 
works by Mcndez & van der Klis (1997) and Kong et al. 
(2000). 

4-2. Estimate of the Innermost Disk Radius 

Assuming that the number of photons emitted from the 
disk is conserved in the Comptonization process, we can 
estimate the innermost radius of the disk from the follow- 
ing equation (Kubota & Makishima 2004), 



photons s ^ cm 



(2) 



^disk ■ 



i^fhc2cosi : 



=0.0165 



(L»/10 kpc) 



Tin Y 

1 kev; 



where -Fjigi^ and F^-^^ are the unabsorbed 0.01-100 keV 
photon flux from the disk and thermal Compton com- 
ponent, respectively. Here it is assumed that the 
Comptonizied emission is isotropic and there are few pho- 
tons that are scattered back to the disk. Table 2 list the 
values of F^-^^y. and F^^^ (for the two Comptonization com- 
ponents) in each epoch; the direct disk flux is found to be 
less than 3% of the total Comptonized one. Assuming 
i = 50° (Section 4.2), we derive the innermost disk radius 
for Epoch-1, -2 and -3 to be ri„ = (98±2)L>8, (107t2^)i:>8, 
(97±2)D8 km, respectively, where Ds is the distance to 
GX 339-4 in the unit of 8 kpc. Here the errors reflect the 
uncertainties in Tin. Multiplying this radius by the cor- 
rection factor, 1.19, for the boundary condition and ratio 
of the color/effective temperature of the multi-color disk 
(see Kubota et al. 1998), yields an actual innermost ra- 
dius of = (117±2)L>8, (128tf )i:>8, (115±2)i:>8 km 
for the three epochs. 
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Fig. 9. Iron K-line fit to the 3-8 keV XIS-O+XIS-3 spectra 
in 2009 March with the diskline model for a fixed emissivity 
index of /3 = —3. The ratio between the data and best-fit 
model are shown in the lower panel. 

4-. 3. Iron-K Line Fitting 

To best constrain the disk geometry from the iron-K 
emission hne profile, we fit the XIS spectra in the 3-8 
keV band with the relativistic disk hne model diskline 
(Fabian et al. 1989) with the reflection continuum self- 
consistently. Here we consider both Ka and K/3 lines; 
although the contribution of the K/3 line is generally weak 
and has been neglected in most of previous work, wc find 
that it affects the final results. Since the parameters can 
be only weakly constrained from the spectra in each epoch 
due to the limited photon statistics, we combined the spec- 
tra of the three epochs separately for XIS-O-l-XIS-3 and 
for XIS-1 to discuss the time-averaged properties over the 
three observations. This is justified because the contin- 
uum fit to the individual broad band spectra gives similar 
(if not the same) results to one another, and the estimated 
inner disk radii (i??°"*) stay constant within 10% level as 
shown in the previous subsection, which is much smaller 
than the statistical error we will obtain from the line pro- 
file analysis. As a continuum, we adopt the results from 
double compPS fit using the broad-band XIS-I-HXD spec- 
tra averaged over the three epochs, and replace the gauss 
model with diskline for the iron line. In the fit of the XIS 
spectra, the normalization and optical depth of the hard 
Comtonization component, the most dominant one in the 
3-8 keV band, are treated as free parameters. Since the 
continuum shape is hard to constrain within the limited 
band coverage, they are allowed to float only within their 
90% errors obtained in the XIS-t-HXD fit, and the other 
continuum parameters are fixed at the best-fit values. 

The diskline model has five parameters, the inner and 
outer disk radius {Rin and -Rout), inclination angle (i), 
power law index of emissivity /3, and intrinsic line center 
energy {Ecen)- We fix i?out = 10000i?g, which is difficult to 
constrain from the data. The line center energies are fixed 
at 6.40 keV (Ka) and 7.06 keV (K/3) to make them con- 
sistent with the low ionization parameter obtained from 
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Fig. 10. Fitting results of the 3-8 keV XIS-O+XIS-3 spec- 
tra with the diskline model, obtained from the observations 
in 2009 March (our data, brighter state) and that in 2008 
September (Tomsick et al. 2009). The top panel plots the 
data (black, with error bars) and the best-fit models of the 
iron-Ka (blue, dashed) and -K/3 lines (blue, dotted) and of 
the continuum (red, solid). The ratio between the data and 
best-fit model are shown in the lower panels for the 2009 and 
2008 data, respectively; the red ones (open circles) correspond 
to the case when the iron-K line is excluded from the model. 
For clarity the XIS-1 data are not plotted. 



Table 3. Parameters of the Iron-K Line Fit* 



Component Parameter 



Value 



diskline (Ka) 


E,,n (keV) 


6.4 (fixed) 




EW (eV) 


87^1? 


diskline (K/3) 


-Beerr (kcV) 


7.06 (fixed) 




Emissivity law j3 


-2.3±0.1 




i (deg) 


46 ±8 




-Rin (-Rg) 


-, o 0+6.4 




intensity ratio (K/3/Ka) 


0.13 (fixed) 



xVd.o.f. 273/271 

Notes. The errors are 90% confidence level for a single parameter. 
* The continuum model is wabs X (diskbb + compPS -f compPS). 
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Fig. 11. Top: Delta chi-squared values plotted against 
obtained from the diskline fit. The dashed, dash-dotted, and 
dotted lines represent 68%, 90%, and 99% confidence limits 
for a single parameter, respectively. Middle: Confidence con- 
tour between the inner radius Rj^ and inclination i at 90% 
(red, outer) and 68% (black, inner) confidence limits for two 
parameters. Lower: Confidence contour between the inner 
radius ijjn and emissivity index /3 at 90% (red, outer) and 
68% (black, inner) confidence limits for two parameters. 



the reflection continuum. We also fix the intensity ratio 
between the K/3 and Ka at 0.13, a theoretical value for the 
iron K-fluorescence process (Kaastra & Mewe 1993). The 
rest are left as free parameters. Because the reflection 
continuum must also be blurred in the same manner as 
the emission line, we tie i?in, Rout, and /3 in the diskline 
model with those of the reflection component included in 
the compPS model. The result on the inclination can be 
compared with the independent constraints presented in 
Figure 1. When the emissivity index is fixed at /3 = — 3, we 
find that the fit is not acceptable (x^/d.o.f. = 323/272), 
showing a significant discrepancy between the data and 
model in both sides of the iron line (see Figure 9). In 
addition, the best-fit inclination angle for /? = — 3 is found 
to be « 20°, which is not reasonable since the black hole 
mass exceeds lOOM© according to Figure 1, unless the in- 
ner accretion disk is heavily warped with respect to the 
binary plane. Thus, we vary f3 in the range — 3 < /3 < — 2 
as a free parameter. 

The result of the fit is plotted in Figure 10, and the 
best-fit parameters are listed in Table 3. We also show 
the increment of in terms of i?in, and confidence con- 
tour between i?in and i, and that between i?in and 13 in 
Figures 11 (a), (b), and (c), respectively. The inclination 
angle now becomes i = 46 ± 8 degrees, which is reasonable 
in terms of the black hole mass (Figure 1). The resulting 
inner radius is i?in = 13.3tg;gi?g, suggesting that the stan- 
dard disk is extended close to the ISCO of a non-rotating 
black hole but is likely truncated at 90% confidence level 
(see Figure 11). 

To compare with i?f°"', the unit of the innermost disk 
radius derived through the iron-K line fit must be con- 
verted to km; i?in = 13.3^6 n-Rg corresponds to i?in = 



200:^go''Mio km, where Mw is the black hole mass in the 



unit of 10 Mq. Thus, the results from the continuum and 
iron-K line profile analysis become consistent if the BH 
mass is the range of (4-16)M0, assuming the distance of 
8 kpc. 

4-. 4- Comparison with 2008 Suzaku data 

Tomsick ct al. (2009) observed GX 339-4 in the fainter 
low/hard state using Suzaku, from 2008 September 24 
22:37:25 to 27 02:19:24 (UT) with a net exposure of 
100 ksec for the XIS and 80 ksec for the HXD (ObsID: 
403067010). The 1-100 keV flux of the source was 
2.4 X 10^^" erg cm^^ sec^^ in this period, which is ap- 
proximately 14 times smaller than that of our 2009 data. 
As described in Section 4.1, the best-fit R^^ value in the 
2009 data is i?in = 13.3^6^0' apparently a factor of >8 
smaller than the 2008 result by Tomsick et al. (2009). We 
note, however, that the derived i?in is coupled with both 
inclination and emissivity index. To eliminate the uncer- 
tainties, we re-analyze the data of the 2008 observation. 
The data reduction is made according to the standard pro- 
cedure from the cleaned event files of processing version 
2.2.11.22. We extract the XIS time- averaged spectra in 
the 0.5-10.0 keV band. The XIS-0 and XIS-3 spectra are 
co-added together with their responses. 

We analyze the spectra with the model composed of an 
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Table 4. Parameters of the Iron-K Line Fit to the Tomsick 
ct al. (2009) Data* 



Component Parameter Value 



Notes. The errors are 90% confidence level for a single parameter. 
* The continuum model is wabs X powerlaw. 

absorbed power law plus two iron line (diskline) compo- 
nents for Ka and K/3 line. For the analysis, we use the 
XIS-O-l-XIS-3 and XIS-1 spectra in the 3-8 keV range. We 
fix the inclination angle at i = 50° and the emissivity index 
at /3 = — 2.3 according to our results of the Suzaku 2009 
data. The quality of the fit is found to be acceptable, 
■y/^ jv = 444/402, with the best- fit photon index of T = 
1.61 ± 0.01 and column density of A^h = 5.31 ± 0.07 x lO^^. 
The XIS spectra, best-fit model, and their ratio are plot- 
ted in Figure 10. Table 4 summarizes the best-fit param- 
eters. It is seen that the iron-K lines in the 2008 data is 
significantly narrower than those observed in 2009. The 
resulting inner disk radius is found to be i?in = 190^90°: 
which is wl4 times larger (at least 5 times larger within 
the statistical errors) than that in our observations. 

5. Near-Infrared Observations and Data Analysis 

We carried out near-infrared photometric observations 
of GX 339-4 on 11 nights over a period from 2009 Feb 27 
to 2009 March 23, using the SIRIUS camera (Nagayama et 
al. 2003) on the 1.4 m IRSF telescope at the South African 
Astronomical Observatory (SAAO). SIRIUS is a JEK^,- 
simultaneous imaging camera, which has three 1024 x 
1024 HAWAII arrays and covers a field of view of 7'. 7 x 7'. 7 
with a scale of 0".45 pixel"^. The log of observations is 
given in Table 5. The typical seeing (FWHM) during the 
observations was 1".4 (3 pixels) in the / band. 

After a standard data reduction procedure for near- 
infrared array images using the IRAF analysis package 
(dark subtraction, flat-field correction and sky subtrac- 
tion), 10 dithered images were combined to reject bad 
pixels and improve the S/N ratio. Hence, each magnitude 
in Figure 12 was determined from 50 or 100 sec integra- 
tion. Aperture photometry was made with an aperture 
of 12-pixel radii. Photometry with smaller apertures gave 
similar results. We used the 2MASS data of about 150 
stars in the field for magnitude calibration. The 1-day av- 
eraged magnitudes in the JHKg band are listed in Table 5. 

The J (1.25 ^m) and (2.14 ^m) light curves with the 
J — Kg color are plotted in Figure 12, together with those 
of Swift/BAT in the 15-50 keV band and RXTE/ASM 
in the 2-12 kcV band. As noticed, the daily averaged 
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Fig. 12. X-ray (15-50 keV and 2-12 keV) and IR (J and Kg 
band) light curves and J — Kb magnitudes of GX 339—4 ob- 
tained with Swift/BAT, RXTE/ASM, and IRSF/SIRIUS in 
2009 March. The abscissa represents the date of MJD-54890, 
where 1 corresponds to 2009 March 1. The IR magnitudes are 
calibrated based on the Two Microns All Sky Survey (2MASS) 
Point Source Catalog, by using stars in the field of view. The 
arrows in each panel show the epochs of our Suzaku observa- 
tions. 

Kg band flux is correlated positively with the hard X-ray 
flux. Significant variation in the J — Kg color by <0.3 
magnitudes is observed on time scale of 100 sec, although 
we find no clear correlation between the color and mag- 
nitude. The H (1.63 fim) band light curve is not plotted 
because they basically follow the same trend as the J and 
Kg light curves. 

The left panel of Figure 13 shows the SED of GX 339- 
4 at the near-infrared to optical wavelengths in the i^Fi, 
form. Utilizing the IRSF data taken between 2009 March 
20 to 24, when the Swift/BAT fiux was nearly constant 
around the Suzaku observation epoch, we plot the aver- 
aged photometry of the J-H-Kg bands, the spectrum with 
the reddest J ~ Kg color, and that with the bluest color. 
For comparison, the near infrared to optical data acquired 
in 1981 and 1997 in the low/hard state are also plotted 
(Corbel & Fender 2002, and references therein). In the 
right panel, a combined SED of GX 339-4 over the radio 
to X-ray bands is shown, although they are not simulta- 
neous except for the IRSF and Suzaku data, which are 
quasi-simultaneous. The interstellar absorption is cor- 
rected in the X-ray spectrum, and all the near infrared 



diskline (Ka) E^^n (keV) 6.4 (fixed) 

EW (eV) 77tll 3 

diskline (K/3) Eccn (keV) 7.06 (fixed) 2 

Emissivity law /3 —2.3 (fixed) | 

i (deg) 50 (fixed) 



Rin (i?g) 190+^^° 



intensity ratio (K/J/Ka) 0.13 (fixed) 
xVd.o.f. 444/402 
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Table 5. The Log of Ncar-infrarcd Observations 



Date* 


Number of Observations 


Integration Time per Frame (sec) 


Averaged Magnitude 








J 


H 


Ks 


Feb 28 


2 


5 


13.43 ±0.03 


12.60 ±0.02 


11.89 ±0.03 


Mar 2 


2 


5 


13.66 ±0.03 


12.80 ±0.02 


12.07 ±0.03 


Mar 3 


4 


5 


13.52 ±0.04 


12.67±0.03 


11.85 ±0.03 


Mar 5 


2 


5 


13.23 ±0.03 


12.40 ±0.02 


11.61 ±0.02 


Mar 7 


9 


5 


13.27 ±0.06 


12.40 ±0.04 


11.63 ±0.05 


Mar 10 


6 


5 


13.21 ±0.05 


12.30 ±0.03 


11.53 ±0.03 


Mar 20 


7 


10 


13.00 ±0.03 


12.12 ±0.02 


11.33±0.02 


Mar 21 


2 


10 


13.03 ±0.02 


12.16±0.01 


11.36 ±0.01 


Mar 22 


6 


10 


12.99 ±0.03 


12.11 ±0.02 


11.38 ±0.02 


Mar 23 


2 


10 


13.05 ±0.02 


12.17±0.01 


11.47±0.01 


Mar 24 


5 


10 


13.16 ±0.03 


12.30 ±0.02 


11.49±0.02 



" All observations were made in 2009. 

t We carried out 10 ditherings at eacli time. 
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Fig. 13. Multiwavclcngth SED of GX 339-4 in the uFu form, corrected for the interstellar absorption/extinction. In the left panel, 
the data from our IRSF/SIRIUS observations, those obtained in 1981 and 1997 (Corbel & Fender 2002) are shown in red, green 
(open triangle), and blue (open circle), respectively. The IRSF fluxes averaged over MJD 54910—54915 are plotted in red, that with 
the reddest J — Ks color are in pink (filled square) , and with the bluest color in orange (open square) . In the right panel, the Suzaku 
spectra and the radio data in 1992-1999 (Fender 2001) and 1997 (Corbel & Fender 2002) arc displayed in red, purple (open square), 
and blue (open circle), respectively, in addition to the same data as the left panel (without the reddest /bluest J — Ks data). The solid 
line shows the estimated contribution of intrinsic multicolor disk emission, including the Compton scattered photons. The dotted 
line represents the upper limit of the contribution from the companion star, assuming a blackbody with an effective temperature of 
4000 K (Zdziarski et al. 2004). The dashed line shows a power law with an index of 1.0 (a = 0.0), fitted to the Fender (2001) results, 
for illustrative purpose. 
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data are deredenned by assuming Ay = 3.7 (Corbel & 
Fender 2002). Both X-ray and infrared fluxes in our ob- 
servations are between those obtained in 1981 and 1997, 
which difi'er by a factor of 4. By fitting the J-H-Kg spectra 
by a power law with oc i^", we obtain a = —0.10 ± 0.08 
(average), which varied between a= — 0.20±0.15 (the red- 
dest data) and a = +0.08 ± 0.10 (bluest). 

6. Discussion 

6.1. Suzaku Results and Inner Disk Structure 

With Suzaku, we have obtained the best quality simul- 
taneous, broad band X-ray spectra ever obtained from 
GX 339-4 in the low/hard state, covering the 0.5-310 keV 
band. We demonstrate that the spectra in the three obser- 
vations can be well described by thermal Comptonization 
of disk photons dominating the flux of the entire band, a 
direct MCD component, and a reflection component with 
an iron-K emission line. The disk photon flux is very small 
(^3% of the Comptonized one), indicating that the inner 
part of the optically thick disk is almost fully covered by 
the corona. The strength of the reflection component (in 
terms of the solid angle, r2/27r « 0.45) and the ionization 
parameter of the reflector (^ ^ 2) are consistent with the 
result from a previous Ginga study in the similar intensity 
state (in 1989 September) analyzed with a power law con- 
tinuum (Ueda et al. 1994), but are now determined with 
much smaller uncertainties thanks to the better energy 
resolution and wider coverage of Suzaku. 

We find that in each epoch, two thermal 
Comptonization components ("double CompPS" model), 
with the same electron temperature (Tc « 175 keV) but 
different optical depths (r « 0.4 and 1). gives a reasonable 
description of the incident continuum spectrum. This 
indicates that the Comptonizing corona can be regarded 
isothermal but has a more complex geometry than that 
assumed in a single zone, spherical model. In fact, 
the corona may have a more disk-like geometry and/or 
inhomogeneous structure, where the optical depths of 
Comptonization as measured from the input source (i.e., 
the standard disk) is not spatially constant. Similar con- 
clusions have been obtained from Cyg X-1 (Makishima et 
al. 2008) and GRO J1655-40 (Takahashi et al. 2008) in 
the low/hard state with Suzaku. Thus, we suggest that 
such physical conditions of the corona may be a common 
feature in the low/hard state of BHBs. Although the 
overall shape of the broad-band spectra of GX 339-4 look 
similar during our observations, the optical depths and 
relative fraction of the two Comptonization components 
are found to be slightly different among the three epochs 
while the electron temperature remains constant at R:!l75 
keV. This suggests that the geometry of the corona varies 
on time scale of several days in the low/hard state even 
at an almost same (< 5%) luminosity level. We note 
that the double CompPS fit still leaves small, systematic 
discrepancies between the data and model in the PIN 
^ 30 keV region. This is not surprising since the change 
of physical parameters of the corona should be continuous 
in reality, while the double CompPS model gives only 
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an approximation. Application of more complex (and 
realistic) Comptonization models is left as future work. 

The analysis of the iron-K line profile has enabled us to 
tightly constrain the inclination angle and innermost ra- 
dius of the optically thick disk. We obtain i ^ 50 degrees, 
which is converted to the black hole mass of 4-16 Mq for 
the allowed range of the mass function of this binary (see 
Section 2). The estimated mass range looks more reason- 
able than huge masses (> 100 Mq) calculated from z = 18 
degrees, if the disk is not strongly warped, as claimed by 
Miller et al. (2006) and Reis et al. (2008) based on their 
detection of a broad iron-K line. The reason for this large 
discrepancy is unclear, but it should be noted that results 
on a very broad iron-K line often strongly depend on the 
continuum modeling and may be partially affected by pile- 
up of the CCD data (Yamada et al. 2009; Done & Diaz 
Trigo 2010). 

The emissivity index obtained from the empirical 
diskline fit also carries important information for under- 
standing the disk structure. In our state, a part of the 
iron-K line emitted from the innermost radii should be 
subject to Comptonization as well, under the situation 
that the hot corona almost fully covers the inner disk. 
This causes significant broadening of the scattered line 
component, with an averaged fractional energy gain of 
AE/E ~ AkTe/mc^ ^ 1, which makes it almost unde- 
tectable as a line. This would partially explain the ob- 
served flatter slope of the emissivity, P = —2.3, than the 
simplest expectation of /3 = —3 when a perfectly plain disk 
is irradiated from a point source with a finite scale height. 
Another reason for the flat slope is that the index will ap- 
proach to /3 = — 2 at outer radii as the scale height of the 
irradiating source becomes negligible compared with that 
of the disk. 

As mentioned in Section 4.3, the innermost radius ob- 
tained from the diskline fit, in units of i?g, are consistent 
with that independently estimated from the MCD contin- 
uum, in units of km, within the uncertainties of the mass 
function and distance. Comparing our i?in value with 
those obtained in previous observations, we can study 
how the inner edge of the accretion disk evolves during 
the low/hard state as a function of luminosity. The 1- 
100 keV flux during our observations is approximately 
14 times brighter than that in the 2008 Suzaku obser- 
vation (Tomsick et al. 2009), which was 2.4 x 10"^° erg 
cm^^ sec^^. As described in Section 4.3, the Rm value 
(i?in = 13.3lg'Q) is at least 5 times smaller within the sta- 
tistical errors than that obtained in our re-analysis of 2008 
data by adopting the common inclination (i = 50°) and 
emissivity index (/3 = —2.3). These results indicate that 
the inner edge moved inward as the luminosity increased 
from 2008 to 2009. 

Table 6 summarizes major previous results on the in- 
nermost radius i?in of the standard disk of GX 339-4 in 
different X-ray luminosities and states, all derived from 
analysis of the iron-K emission line with a relativistic disk- 
line model. We also list other line parameters including 
the inclination and emissivity index, with which R^ could 
be coupled. The original results obtained by Miller et al. 
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Table 6. Summary of Previous Results of Iron-K Line Fit Using a Relativistic Disk-line Model. 
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(2006) in the low/hard state and Miller et al. (2008) in the 
very high state have been recently re-examined by Done & 
Diaz Trigo (2010) and Yamada et al. (2009), respectively, 
and it is pointed out that both data were significantly 
affected by pile-up. The re-analysis show that the data 
favor larger values of Rm in both cases, although the new 
constraints become weaker due to poorer photon statis- 
tics. 

Figure 14 plots Rin with its 90% statistical error for 
a single parameter, as a function of Eddingtion fraction 
Lx/^Edd {D = 8 kpc and black hole mass of IOMq are 
assumed). This figure indicates that the standard accre- 
tion disk of GX 339-4 evolves inward as the luminos- 
ity increases, with a significant change in the range of 
^ 0.001 < ix/isdd 0.02 when the source remains in 
the low/hard state, as concluded by Tomsick et al. (2009). 
Our result in the bright low/hard state in 2009 is consis- 
tent with that of Done & Diaz Trigo (2010) obtained at 
a similar luminosity if R^ ~ 20i?g. This suggests that at 
ix/^Edd = 0.02 the standard disk is likely truncated be- 
fore reaching to 6i?g, the ISCO for a non-spinning black 
hole, although our Suzaku result alone cannot rule out the 
possibility that the inner edge lies below 6i?g at 99% confi- 
dence. If GX 339-4 has an extremely-rotating black hole 
as claimed by Miller et al. (2008), our 2009 result then 
would be strong evidence that the disk docs not extend to 
the ISCO in the bright low/hard state. Further observa- 
tions with similar or better energy resolution not limited 
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Fig. 14. Estimated innermost radii of the standard disk 
of GX 339-4 derived from our results from the 2008 
(left) and 2009 (right) Suzaku data with an inclination 
of i = 50° and an emissivity index of /3 = —2.3, plot- 
ted against fraction of Eddington luminosity (ix/^Edd)- 



by the photon statistics and systematic errors caused by 
pile-up would be necessary to further clarify the right loca- 
tion of the inner edge of the disk in the bright (■~0.01iEdd) 
low/hard state. 
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6.2. IRSF Results and Jet Energetics 

We examine the origin of the multiwavelength SED 
of GX 339-4 in the low/hard state with different emis- 
sion components, utihzing our quasi-simultaneous near in- 
frared and X-ray spectra combined with previous data. In 
the right panel of Figure 13, we plot the estimated contri- 
bution of the "intrinsic" MCD component including pho- 
tons that are Compton-scattered in the corona, through 
the photon number conservation, assuming that the disk 
is extending to infinity (for illustrative purpose) with the 
temperature proportional to r~^/^, where r is the radius 
from the black hole. The results are based on the double 
CompPS fit, which gives the innermost disk temperature 
of 0.22 keV. We also plot the upper limit of the contribu- 
tion from the companion star, assuming a blackbody of 
4000 K. Here we refer to the r-band magnitude of > 21.4 
by Zdziarski et al. (2004), which is then corrected for ex- 
tinction corresponding to Ay = 3.7. As noticed from the 
figure, the blackbody component of the multi-color disk 
and companion star contributes only < 0.5% and ^ 1% 
of the total fiux in the near-infrared range, respectively. 
In the radio band, a slightly inverted power law spectrum 
is observed with an energy index of a «0.1-0.2 (Corbel 
et al. 2000), which is interpreted as an optically thick 
synchrotron emission from the compact jets, a common 
feature observed from BHBs in the low/hard state. 

Our averaged IRSF spectrum shows a somewhat smaller 
energy index, a ~ —0.1, than the radio spectrum; the cor- 
responding J — Ks color is similar to that of the 1981 ob- 
servations compiled by Corbel &: Fender (2002). As shown 
by them, this color is consistent with a superposition of 
two different components with a = —0.6 and a = 2.1, which 
can be explained by the optically thin synchrotron radia- 
tion from the jets and the reprocessed, thermal emission 
by X-ray irradiation from outer parts of the disk, respec- 
tively (see also Markoff et al. 2003; Gandhi et al. 2010). It 
is thus indicated that the transition point from optically 
thick to thin regime of the synchrotron emission, defined 
as Ut, must appear at frequencies just below or around the 
near-infrared band, ~ lO^"' Hz. Note that this conclusion 
does not match the SSC model discussed by Coriat et al. 
(2009) for the low/hard state of GX 339-4 where the H 
band corresponds to an optically thick emission in our lu- 
minosity range [L / ^Edd ~ 0.01 for A/bh = lOM© and = 8 
kpc). The rapid variability observed in the optical band 
(e.g., Gandhi et al. 2008) suggests that the power law com- 
ponent of synchrotron radiation extends there, although it 
is not clear at which frequency the spectrum breaks due 
to radiative cooling. In our picture, it should be below 
the X-ray band, since the X-ray fiux must be dominated 
by thermal Comptonization that can well account for the 
Suzaku spectra. 

The observed synchrotron luminosity and the limit for 
i^t can be used to constrain the physical parameters of 
the jet base, from which most of the observed near in- 
frared luminosity is emitted according to the compact jet 
model (Blandford & Kdnigl 1979). For simplicity, here 
we consider a single zone model representing this region. 



We neglect any effects of relativistic beaming from the 
bulk motion of the jet, since the Doppler factor is esti- 
mated to be 5 pal. 03 for an inclination of 50° and an as- 
sumed intrinsic jet velocity of 0.9c. Let us assume that 
the electron number density at the Lorentz factor 7 is 
given as n{'~f) = A'y~P (for p > 2) in the range between 
7min < 7 < 7max- The syuchrotrou luminosity density in 
an optical thin region is expressed as (see e.g., Rybicki & 
Lightman 1979) 

" 2tt{p+1) V4 12/ U 12/ 

3-p 

/ ly \ ^ P+i 
X VAuB(TTc[ — (sin0) ^ ^ 

where V is the volume {V = AttR^ /3 for a sphere with a 
radius of R), ub = B'^ /{Sir) is the energy density of the 
magnetic field, ctt is the cross section of Thomson scat- 
tering, vb = eB / {2'Kmc) is the gyro frequency (e and m 
are the electronic charge and the electron mass, respec- 
tively), r(y) is the gamma function of argument y, and 9 is 
the pitch angle, which is the angle between magnetic field 
and electron velocity. By assuming equipartition of energy 
density between the electrons' kinetic energy and the mag- 
netic field, that is A= (p — 2)ub/ {mc^) with 7inin = 1, we 
can estimate the strength of the magnetic field from an 
observed luminosity for given R and p. Further, using an 
strict form of the synchrotron absorption coefficient, 

(p-2) ^P+3 f3p + 2\ /3p + 22\ 

2 O-T -3 I V \ 2 , . ^. E+l 

Uj ' ' 

the frequency vt below which the synchrotron emission 
becomes optically thick for a region size of R can be deter- 
mined from the condition a^i? = 1. Taking i^L^, « 3 x 10'^^ 
erg s-^ at i/= 1.4 x 10" Hz (for D^8 kpc), p = 2.2 (corre- 
sponding to the spectral index a = (1 — p)/2 = —0.6), and 
i^t ^ 10" Hz, we estimate B « 5 x lO'' G and i? « 6 x 10* 
cm, which corresponds to ^ 4 x 10^i?g. More precisely, 
dependence of B on R and that of i/t on B and R are 
given as 

^-5x10^ ( ( \ G 

V3 X 1035 erg s-iy \6 x 10^ cm J 

and 

p+e 2 

j^t « 1 X 10" ^-r- ^^5 Hz, 

VSxlO'iG/ Vex 108 cm/ 

respectively, where the terms of the pitch angle are aver- 
aged (i.e., integrated over solid angle and normalized by 
An). For this B value, the corresponding energy density of 
the magnetic field of 8.4 x 10'' erg cm"^ and the Lorentz 
factor of electrons emitting the v — l.Ax 10^^ Hz photons 
is 7 « 40. Our results hold within a factor of 2 when wc as- 
sume "particle-dominant" energy balance with ~ lOw^, 
as observed from jets of blazars (e.g., Inoue & Takahara 
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1996). The magnetic-field strength in the jet base is con- 
sistent with the Giannios (2005) result ol _B = 10^ G ob- 
tained for XTE J1118-t-480 in the low/hard state from 
modelling of the multi- wavelengths SED. 

Based on these parameters, we roughly evaluate the ex- 
pected contribution of Compton scattering by the non- 
thermal electrons, which could produce X-ray and 7-ray 
photons from seed photons of the synchrotron radiation 
(synchrotron self Compton; SSC) and external photons 
(external Compton; EC) from the accretion disk, respec- 
tively, for a typical Lorentz factor of 7 = 40. The energy 
density of synchrotron photons in the near infrared to op- 
tical bands is estimated as 

-"sync ^ isync/(47ri?^c) « 7 X 10^ crg cm"^, 

where igync, which depends on 7max, is assumed to be 
Iq36 gj.g g-i_ Siniilarly, that of external photons observed 
in the X-ray to Gamma-ray bands is 

wcomp ^ icomp/(47r/^c) « 2 X 10^ erg cm"^, 
where Lcomp ~ 2 x 10^^ erg s~^ and the distance from the 
disk and Comptonized corona to the jet base is assumed 
to be / = 5 X 10^ cm that corresponds to the time delay 
(w 150 ms) between the optical and X-ray cross correlation 
signals (Gandhi et al. 2008). 

Thus, the photon energy density of these seed photons is 
wlO-40 times smaller than that of the magnetic field, ub~ 
8.4 X 10^ erg cm"''. As noticed from Figure 13, recalling 
that the ratio of the luminosity between synchrotron and 
SSC by the same electrons is proportional to us/wsync, 
we estimate that the contribution of Comptonized pho- 
tons observed in the 1-100 keV X-ray band (SSC) should 
be approximately 0.4% of the total observed luminosity, 
which is dominated by thermal Comptonization by the hot 
corona. This justifies our X-ray spectral model, which ig- 
nores any non-thermal emission from the jets. The EC 
component, however, may contribute to the high energy 
emission above '-^l MeV. 

Similarly, we also ignore the jet synchrotron emission as 
seed photons for thermal Comptonization by the corona 
in our model, because the emitting regions in the jets are 
assumed to be distant from the corona and disk. However, 
we cannot rule out models where the corona and jet base 
are co-spatial and a Comptonized component of the syn- 
chrotron photons significantly contributes to the X-ray 
emission, as discussed in Markoff et al. (2005). We leave 
it for future work to apply self-consistent jet and corona 
models to the simultaneous, multi-wavelengths SED in- 
cluding Suzaku data. 

To examine the correlation between the Kg and hard 
X-ray light curves (Figure 12), we plot the one-day aver- 
aged 15-50 keV flux (Fx) versus Kg band flux {Fks) in 
Figure 15. The data point taken from the 1981 dataset 
(Corbel & Fender 2002) is also plotted. By fitting IRSF 
and Swift points with a power law of Fks oc F£ , we ob- 
tain F = 0.45 ± 0.06. This value is close to the Coriat 
et al. (2009) resuh from GX 339-4 (0.47 < F < 0.49), 
which was obtained between the H and 3-9 kcV bands. 
Similar correlation is reported between the radio and X- 
ray fluxes with a steeper slope of 0.7-0.9 in the low/hard 




10-9 2x10-9 5x10-9 10-* 

X-ray 10-50 keV flux (erg/cmVsec) 

Fig. 15. The relation between the X-ray (15—50 keV) and IR 
(Ks) fluxes of GX 339-4 with Swift/BAT and IRSF/SIRIUS 
in 2009 March. The result of the 1997 observation (Corbel 
& Fender 2002) is also included (blue, open square). Each 
Kg band flux is calculated by averaging all the magnitudes 
obtained on the same day. The data points except for the 
1997 data are fitted by a power law with Fks <x F^)- The 
best-fit result (7^ = 0.45) is shown as a solid line in this figure. 

state(Corbel 2003). These correlations confirm the strong 
link between the jet formation and accretion flow, as di- 
rectly indicated by the strong cross correlation between 
the optical and X-ray bands on much shorter (< 1 sec) 
time scale for this source. 

The slope of the correlation and its dependence on the 
wavelength and luminosity give us important clues to un- 
derstand the origin of the SED. As mentioned above, our 
result is not compatible with models where the near in- 
frared emission is assumed to be an optically thick syn- 
chrotron emission. We note that, however, the thermal 
emission from the irradiated disk must be taken into ac- 
count as well. Its spectra is represented as a superposition 
of blackbody radiation that results in a flat SED if the 
conversion efficiency from the irradiated flux is constant 
over the whole disk (Frank et al. 2002). The relation be- 
tween the incident luminosity (Lx) and that of thermal 
emission from an irradiated outer disk (-Ljy.irr) is given by 
ii/,iir oc Lx at frequencies corresponding to the flat part, 

and ijy.Rj oc in the Rayleigh- Jeans region (see Coriat 
et al. 2009). Thus, it is possible that a contribution from 
the irradiate disk in the Rayleigh- Jeans part could make 
the observed slope F somewhat flatter than those expected 
from jet models. To prove the origins of the near IR emis- 
sion, detailed modeling of the complete SED would be 
necessary, which we leave for our future investigation. 

7. Conclusion 

We observed GX 339-4 with Suzaku and IRSF in the 
low/hard state at 2% of the Eddingtion luminosity. The 
conclusions are summarized as follows. 

1. The broad-band Suzaku spectra are well represented 
by a thermal Comptonization by hot ~ 200 keV elec- 
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trons of seed photons from the optically thick disk, 
with a small contribution of the direct disk compo- 
nent. The corona has at least two different optical 
depths, r w 0.4 and 1, indicating inhomogeneous 
structure. 

2. Analysis of the iron-K line with a diskline model 
yields an innermost radius of i?in = (13.3tg;g)i?g 
with an estimated inclination of ~ 50° . 

3. The Suzaku results indicate that the optically thick 
accretion disk is truncated before the ISCO and its 
inner regions are almost fully covered by the hot 
corona. The inferred innermost radii from the con- 
tinuum fit and the iron-K line profile are consistent 
each other. 

4. The near infrared fluxes are correlated with that of 
hard X-rays by the relation Fks oc F^ "*^, with vari- 
able spectral indices from -0.20 to -1-0.08. The spec- 
trum can be explained by a sum of an optically thin 
synchrotron emission and a thermal emission from 
the X-ray irradiated outer disk. We estimate the 
magnetic field and size of the jet base to be 5 x 10'* 
G and 6 x 10* cm, respectively. 
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